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A.  STATEMENT  OF  THE  PROBLEM  STUDIED 

This  research  develops  a  new  analytical  and  numerical  method,  namely  the 
Distributed  Transfer  Function  Method  (DTFM),  for  the  modeling  and  analysis  of 
general  complex  flexible  systems  that  are  assemblages  of  multiple  continua  and 
lumped  parameter  systems.  The  research  introduces  a  distributed  transfer  function 
formulation,  based  on  which,  the  eigensolutions,  static  and  dynamic  response  of 
complex  flexible  systems  subject  to  a  variety  of  forcing  sources,  and  damping, 
inertial  and  gyroscopic  effects  are  thoroughly  investigated.  The  distributed  transfer 
function  method  is  also  used  to  study  control-structure  interactions  which  are 
essential  to  active/passive  structural  control  and  smart  structures.  The  conducted 
research  is  fundamental  to  the  dynamics  of  structures  and  machines,  and  has 
important  applications  in  optimal  design  of  military  structures  and  equipment. 
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B.  SUMMARY  OF  RESEARCH  RESULTS 

The  objective  of  this  research  is  to  develop  an  innovative  tool,  namely  the  Distributed  Transfer 
Function  Method  (DTFM),  for  modeling,  analysis  and  optimal  design  of  complex  flexible  militaiy 
structures.  The  research  project  has  gone  through  the  four  stages  shown  in  Fig.  1:  fundamental 
development  of  DTFM-based  synthesis;  modeling  of  flexible  components;  static  and  dynamic 
analysis  of  complex  structural  systems,  and  study  on  control-structure  interaction.  In  sequel,  the 
results  obtained  from  each  stage  will  be  detailed. 


Fig.  1 
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1.  Fundamental  Development  of  DTFM-Based  Synthesis 

The  DTFM-based  synthesis  developed  follows  four  steps  as  illustrated  in  Fig.  2.  The  transfer 
function  synthesis  deals  with  complex  flexible  systems  (CFS)  with  multiple  branches  and  closed- 
loops,  distributed  connection  of  flexible  bodies,  and  multi-point  connection  of  flexible  and  lumped 
subsystems.  The  synthesis  combines  the  advantages  of  both  the  transfer  function  method  and  the 
finite  element  method,  namely  the  exact  and  closed  form  of  the  distributed  transfer  function 
method,  and  the  flexibility  of  the  finite  element  method  in  dealing  with  complex  flexible 
subsystems  [Jl]^.  Theoretical  and  numerical  work  reveals  the  following  important  properties  of 
the  proposed  synthesis: 

(a)  The  synthesis  provides  exact  and  closed-form  solutions  for  many  CFS,  without  the  need 
for  discretization. 

(b)  The  synthesis  does  not  need  to  select  particular  shape  functions  for  specific  system 
parameters  and  boundary  conditions. 

(c)  The  synthesis  is  efficient  in  treating  various  constraints,  and  multi-body  coupling. 

(d)  The  synthesis  can  precisely  predict  high  frequency  responses  of  CFS. 

(e)  The  synthesis  is  highly  accurate  in  estimating  high-gradient  stress  distributions  of  CFS. 


Step  1.  Decompose  the  structures  into  a 
number  of  basic  subsystems 

^  iAvA 

Step  2.  Obtain  the  analytical/semi-analytical 
and  closed-form  transfer  functions  for  each 
subsystem 

i 

Distributed  T  ransfer  Functions 

/  \ 

4  4  4 

1  1  T 

Resportse  External  Boundary 

Disturbance  Disturbance 

Step  3.  Assemble  the  subsystems  => 

Dynamic  equilibrium  equation 

K{s)uis)  =  q{s) 

Step  4.  Solve  for  u{s) 

static  and  dynamic  response 
mechanical  energy  and  power  flow 
eigensolutions,  stress  distributions 

Fig.  2 

One  unique  feature  of  the  distributed  transfer  function  formulation  is  that  different  flexible 
systems  are  systematically  treated  with  the  same  algorithm.  Inertia,  stiffness,  axial  loads, 

^  [J1]  stands  for  the  reference  J1  listed  in  the  publication  section  in  pp.  11-12  of  this  report 
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gyroscopic  and  centrifugal  forces,  circulatory  effects,  and  various  boundary  and  external 
disturbances  are  conveniently  described  without  requiring  different  derivation,  which  is  often  the 
case  with  other  modeling  techniques.  This  feature  allows  computerized  symbolic  manipulation  for 
general  multi-body  systems,  and  makes  the  proposed  modeling  method  flexible  in  dealing  with  real 
engineering  problems. 

Computer  programs  for  implementing  the  synthesis  have  developed.  These  programs  are 
useful  in  study  of  the  eigenvalue,  response,  stability  and  control  problems  of  CFS. 

2.  Modeling  of  Complex  Flexible  Systems 

One-dimensional  flexible  subsystems  [J2] 

Determination  of  the  distributed  transfer  functions  of  flexible  subsystems  is  essential  to  the 
proposed  transfer  function  synthesis.  The  research  first  obtains  exact  closed-form  distributed 
transfer  functions  of  uniform  one-dimensional  flexible  subsystems  such  as  strings,  bars,  beams, 
and  beam-columns,  and  then  extends  the  derivation  to  non-uniform  systems  such  as  beams  with 
varying  cross-section  and  mass  density.  Two  numerical  methods  are  developed  to  estimate  the 
transfer  functions  of  non-uniform  systems.  Compared  to  the  existing  methods  the  transfer  function 
modeling  is  more  accurate  and  efficient. 

Constrained/combined  flexible  systems  [J3,  T2]^ 

The  research  investigates  flexible  systems  that  are  constrained  by  springs  and  dampers,  and  are 
combined  with  rigid  bodies.  This  class  of  systems  has  wide  applications  in  army  equipment  and 
facilities  such  as  robotized  armament  systems,  rotorcraft,  missiles  and  ground  vehicles.  A 
systematic  modeling  and  analysis  method  for  constrained/combined  flexible  systems  has  been 
developed.  In  the  method,  constraint  forces  and  flexible-rigid  body  coupling  are  analyzed  using  the 
obtained  distributed  transfer  functions  and  transfer  matrices  for  the  constraints  and  rigid  bodies. 
Formulas  for  calculating  eigensolutions  and  frequency  response  are  derived.  One  advantage  of  the 
developed  method  is  that  constraints  and  combined  lumped  systems  are  specified  in  a  simple  and 
compact  matrix  form  which  is  convenient  for  computer  coding  and  numerical  simulation. 

As  application  of  the  transfer  function  method,  four  types  of  systems  have  been  considered:  (a) 
coupled  string-rigid  body  systems;  (b)  beams  connected  with  rigid  blocks;  (c)  beam  frames  and 
truss  systems;  and  (d)  a  multi-body  vehicle  structure.  The  main  concern  is  on  the  coupling  and 
interaction  of  flexible  and  rigid  bodies.  Computer  programs  for  calculating  natural  frequencies  and 
mode  shapes,  and  frequency  response  are  being  developed  using  MATLAB.  Numerical  analysis 
and  comparison  of  the  proposed  method  with  the  existing  modeling  techniques  show  many 
advantages  of  the  transfer  function  method.  It  is  believed  that  the  computer  programs  developed  in 
this  research  will  be  quite  useful  for  real  army  applications. 

Complex  flexible  rotor  systems  [J16,  Tl] 

The  distributed  transfer  function  synthesis  is  also  applied  to  rotor  systems  assembled  from 
multiple  flexible  and  rigid  components,  and  constrained  by  bearings  and  gears;  see  Fig.  3.  The 
method  delivers  highly  accurate  and  closed-form  analytical  solutions,  and  is  capable  of  treating 

^  [T2]  stands  for  the  reference  T2  listed  in  the  publication  section  in  pp.  11-12  of  this  report 
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non-self-adjoint  effects  of  gyroscopic  and  damping  forces,  general  boundary  conditions  and  multi¬ 
body  coupling.  It  is  shown  that  the  proposed  method  provides  a  useful  analysis  tool  for  many 
problems  in  rotor  dynamics. 


Co-axial  rotors 


5  6  8  9  10 

Fig.  3 


Stepped,  ring-stiffened  and  laminated  cylindrical  shells  [J4,  J6,  J8] 


Stepped  and  ring-stiffened  cylindrical 
shells,  and  thick  and  laminated  shells  have 
been  modeled  analyzed,  with  interest  in  tank 
cannons  and  other  weapon  systems.  The 
DTFM  is  combined  with  Fourier  series  ex¬ 
pansion  to  yield  highly  accurate  pre-dictions 
of  shell  response,  and  dis-continuities  in  the 
distributions  of  shell  internal  forces.  An 
example  is  shown  in  Fig.  4  where  the 
internal  forces  of  a  cylindrical  shell  of  two 
stiffeners  in  free  vibration  are  plotted  along 
the  shell  longitudinal  direction. 


X 

Fig.  4 


Arbitrarily  shaped  plates  and  2-D  elastic  regions  [J7,  JIO,  J15] 

Extending  from  analysis  of  one-dimensional  continua,  the  proposed  DTFM  has  found  useful 
and  important  applications  in  modeling  of  complex-shaped  plates  and  2-D  elastic  regions.  Unlike 
the  cylindrical  shells,  the  complex  geometry  and  boundary  conditions  of  these  2-D  continua  do  not 
permit  Fourier  series  expansion.  To  overcome  the  difficulty,  a  generalized  technique  called  the 
Strip  Distributed  Transfer  Function  Method  (SD'l'FM)  is  developed  for  the  two-dimensional 
continua.  In  this  semi-analytical  method,  a  multi-body  region  is  first  decomposed  into  a  number  of 
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basic  2-D  subregions  (rectangular  or  sectorial).  Each  subregion  is  then  divided  into  a  number  of 
strips;  the  response  of  each  strip  is  interpolated  in  the  unknown  nodal  line  displacements,  which 
are  functions  of  the  strip  longitudinal  coordinate  and  time.  The  nodal  line  displacements  are 
determined  in  an  exact  and  closed  form  by  the  distributed  transfer  functions  that  are  defined  along 
the  strips.  Synthesis  of  the  subregions  using  the  strip  distributed  transfer  functions  yields  semi- 
exact  and  semi-analytical  prediction  of  the  static  and  dynamic  response,  natural  frequencies  and 
buckling  loads  of  the  multi-body  continuum. 


The  SDTFM  is  shown  to  be  much  more  accurate  and  efficient  than  the  Finite  Element  Method 
(FEM).  For  instance,  consider  the  calculation  of  the  natural  frequencies  of  a  cantilever  plate  in 
transverse  vibration  (Fig.  5a),  and  the  stress  distribution  of  an  L-shaped  2-D  elastic  region 


a/2 

(b) 

— 

P 

b 

B 

A 

Fig.  5 


(Fig.  5b),  on  line  A-B  where  the  stress  <Sy  becomes 
singular  at  point  A.  Table  1  and  Fig.  6  show  the  high 
accuracy  of  SDTFM  in  predicting  higher-mode 
eigenvalues  and  high-gradient  stresses:  with  just  4 
and  12  strips  in  the  two  cases,  respectively,  SDTFM 
delivers  the  same  accurate  results  as  obtained  by 
FEM  with  900  and  2,700  elements. 


Table  1.  N(»-dimensional  natural  frequencies  of  a 
cantilever  plate  (a/b  =  1.5) 


Mode 

No. 

SDTFM 

4  strips 

FEM 

4x4 

FEM 

8x8 

FEM 

30x30 

1 

1.5352 

1.5314 

1.5346 

1.5348 

2 

5.1825 

5.2300 

5.1970 

5.1808 

3 

9.5422 

10.252 

9.7418 

9.5534 

4 

17.489 

18.573 

17.854 

17.502 

5 

23.820 

25.081 

24.213 

23.812 

10 

58.285 

64.454 

60.252 

58.363 

15 

87.625 

94.502 

91.959 

88.364 

20 

124.80 

157.38 

137.71 

125.84 

The  SDTFM  is  not  limited  to  rectangular  and  sectorial  shapes, 
it  is  valid  for  elastic  region  with  arbitrarily  curved  boundaries. 
Through  a  coordinate  transformation,  an  arbitrarily  shaped  spatial 
domain  of  the  continuum  is  mapped  onto  a  rectangular  region, 
where  the  isoparametric  strip  distributed  transfer  functions  of  the 
plate  are  introduced.  Consequently,  semi-exact  and  semi-analytical 
solutions  for  various  static  and  dynamic  problems  can  be 
systematically  obtained.  As  an  example,  consider  the  plate  in  Fig. 
7,  where  the  two  curved  boundaries  are  described  by 
flo(y)  =  -0-lsin37ty  and  ai(y)  =  -1.2y(y-l)-l-0.2.  The 
calculated  nondimensional  natural  frequencies  listed  in  Table  2 


Fig.  7 
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indicates  that  only  a  few  strips  are  needed  in  SDTFM  to  obtain  the  results  that  require  thousands  of 
elements  in  a  finite  element  analysis.  Analysis  of  other  shaped  plates  has  been  considered,  and  also 
reveals  the  high  accuracy  and  efficiency  of  the  method. 


Table  2.  Natural  frequencies  of  the  plate  with  two  different  curved  boundaries 

(ns  -  number  of  strips) 


Method 

Xi 

^2 

X13 

X4 

>^5 

SDTFM 

ns  =  2 

3.0008 

5.5244 

11.3428 

15.8556 

16.4639 

ns  =  4 

2.5970 

4.9815 

9.1872 

14.9220 

14.9308 

ns  =  6 

2.5428 

4.7922 

8.6906 

13.9231 

14.7647 

ns  =  8 

2.5326 

4J326 

8.5532 

13.5813 

14.7411 

FEM  (72x  50  elements) 

2.5407 

4.8422 

8.6978 

13.3843 

14.8493 

3-D  elastic  bodies  [J17] 

The  concept  of  the  strip  distributed  transfer  function  is  further  extended  to  three-dimensional 
elastic  solids.  Semi-analytical  solutions  are  obtained  using  prism  distributed  transfer  functions, 
which  are  defined  on  prisms  of  the  elastic  region,  instead  of  strips.  Again,  one  unique  feature  of 
the  method  is  to  be  able  to  model  multi-body  continua  of  arbitrary  boundaries,  using  a  smaller 
number  of  prisms. 


Visco-elastic  systems  [J18] 

The  DIPM  developed  is  extremely  convenient  for  modeling  and  analysis  of  structures  made  of 
viscoelastic  material,  or  with  viscoelastic  damping  layers.  For  a  linear  viscoelastic  continuum  (1-D, 
2-D  or  3-D)  its  constitutive  relation  in  the  Laplace-transform  domain  can  be  described  by 
o(x,j)  =  sR(x,s)tix,s),  where  O,  R  and  e  are  the  Laplace  transforms  of  the  stress,  strain  and 


relaxation  function  of  the  continuum,  and  s  is  the 
Laplace  transform  parameter.  By  making  use  of 
the  above  .r-domain  constitutive  relation,  DTFM 
can  directly  model  and  analyze  the  viscoelastic 
structure  without  the  need  for  extra  effort.  As  a 
result,  the  proposed  method  avoids  dealing  the 
time-domain  convolution  integrals  that  are 
commonly  encountered  in  many  existing  analysis 
technique.  Furthermore,  through  a  specific  inverse 
Laplace  transform  procedure,  closed-form 
transient  response  of  a  viscoelastic  structure  is 
obtained  for  the  first  time. 


Besides  viscoelastic  materials  or  damping  layers,  the  DTFM  synthesis  is  also  capable  of 
modeling  structures  with  piezoelectric  layers.  Indeed,  the  method  developed  indeed  provides  a 
useful  tool  for  modeling  and  analysis  of  smart  structures  hke  the  hypothetic  one  in  Fig.  8. 
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3.  Static  and  Dynamic  Analysis 

As  shown  in  Fig.  1,  the  SDTFM-based  synthesis  of  a  complex  flexible  system  leads  to  the 
dynamic  equilibrium  equation 

K(s)u(s)  =  q{s)  (1) 

where  u(s)  is  the  vector  of  the  displacements  at  the  points  (called  nodes)  where  flexible 
subsystems  are  interconnected,  q(s)  is  the  vector  of  effective  loads  at  the  nodes,  and  K{s)  is  a 
transcendental  dynamic  stiffness  matrix  which  is  assembled  from  the  distributed  transfer  functions 
of  all  the  subsystems.  Based  on  the  above  equation,  analytical  or  semi-analytical  solutions  of 
various  static  and  dynamic  problems  have  been  obtained. 

Eigenvalue  solutions  [Jl,  J2,  J4,  JIO,  J15] 

For  the  flexible  system  in  free  vibration,  its  characteristic  equation  is  detK(s)  =  0  from  which 
the  system  eigenvalues  k=  1,  2,...,  can  be  determined.  Once  the  eigenvalues  are  known,  the 
corresponding  eigenfunctions  are  obtained  in  two  steps:  first,  determine  the  non-trivial  solution 
W  of  K(X,;fc){<|)jt}  =  0;  second,  substitute  {<t>fc}  into  the  transfer  function  representation  of  the 
response  of  each  subsystem  to  give  the  mode  shape  distributions. 

In  determining  the  eigenvalues  (eigenfrequencies  and  damping  ratios,  or  buckling  loads)  of  the 
flexible  system,  one  needs  to  solve  the  transcendental  characteristic  equation  detK(s)  =  0.  A 
numerical  root  searching  scheme  has  been  developed  based  on  explicit  estimation  of  the  derivative 
d{detK(s)}/d[s.  In  addition,  to  assure  fast  convergency  of  the  scheme,  good  initial  guesses  of  the 
roots  are  made  by  using  the  eigenvalue  inclusion  principles  and  the  spectrum  shift  approach. 

Static  and  frequency  response  [Jl,  J7,  J8,  JIO] 

In  static  analysis,  the  Laplace  transform  parameter  s  in  Eq.  (1)  is  set  zero.  For  the  flexible 
system  under  a  harmonic  excitation  of  frequency  ca,  its  dynamic  response  is  obtained  by  setting 
s  =  i(D,  I  =  V^.  in  Eq.  (1).  Thus,  closed-form  estimation  of  static  and  dynamic  displacements 
and  stresses  is  systematic  and  efficient,  as  has  been  shown  in  Fig.  6. 

Transient  response  analysis  [J3,  J9,  Jll,  J12,  J18] 

An  irmovative  method  for  evaluation  of  transient  response  of  a  complex  flexible  system  with 
damping  and  gyroscopic  effects,  of  viscoelastic  and  smart  materials,  and  under  active  control 
actions  has  been  developed.  This  method  obtains  closed-form  inverse  Laplace  transform  of  the  s- 
domain  response  by  using  the  modal  information.  Two  unique  features  of  the  method  are:  (i)  the 
method  does  not  need  orthogonality  relations  for  the  system  eigensolutions,  which  actually  are 
difficult  to  find  for  the  flexible  systems  in  consideration;  and  (ii)  the  method  does  not  assume  the 
completeness  of  the  eigenfunctions  in  a  function  space. 

Non-self-adjoint  effects  [J2,  J5,  J9,  J16,  C6,  C7]^ 

In  case  study,  the  research  investigated  the  effects  of  non-proportional  damping,  Coriolis 
acceleration,  impact,  follower  forces,  and  inertial  forces  on  the  dynamic  behaviors  of  complex 
flexile  systems.  The  stability  of  the  flexible  system  under  those  effects  are  predicted  through  use 

^  [C6]  stands  for  the  reference  C6  listed  in  the  publication  section  in  pp.  11-12  of  this  report 
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of  the  closed-form  frequency  response.  The  convenience  and  accuracy  of  using  the  distributed 
transfer  functions  is  obviously  observed. 


4.  Control-Structure  Interaction 


As  have  shown  in  the  previous  sections,  the  DTFM  is  accurate  and  efficient  in  modeling  and 
analysis  of  complex  flexible  systems.  The  focus  of  this  part  of  research,  therefore,  is  on  the 
usefulness  of  the  DTFM-based  synthesis  in  design  of  active/passive  controllers  for  complex 
flexible  systems,  and  in  design  of  smart  structures.  Full  development  of  sophisticated  control 
algorithms  is  beyond  the  scope  of  this  project.  However,  to  serve  as  a  precursor  to  future 
investigations,  the  following  three  issues  have  been  selected  in  the  case  study. 

4.1.  Control  system  formulation  by  distributed  transfer  functions  [Jl,  J12,  C6,  Tl] 

The  flexible  system  in  consideration  consists  of  multiple  distributed  subsystems  (e.g.,  flexible 
components,  damping  treatment,  and  smart  material  layers)  and  lumped  subsystems  (e.g.,  rigid 
bodies,  oscillators,  dampers,  and  finite-dimensional  feedback  controllers).  Mathematically,  the 
system  response  is  described  by  a  set  of  coupled  partial  and  ordinary  differential  equations  with 
complicated  boundary  and  matching  conditions.  Conventional  discretization  techniques  like  FEM 
usually  leads  to  a  mo^l  of  a  large  number  of  degrees  of  freedom,  which  usually  has  computational 
difficulties  in  on-line  control.  In  this  research,  it  is  shown  that  for  a  given  frequency  spectrum,  an 
accurate  model  of  a  few  modal  coordinates  can  be  easily  obtained  from  the  distributed  transfer 
function  formulation.  With  that  model,  simple  and  effective  controllers  can  be  convenient 
designed.  Furthermore,  the  study  indicates  that  the  DTF-based  control  system  formulation  for 
high-frequency  dynamics  does  not  require  increased  number  of  unknowns  or  computation  time, 
which  is  an  attractive  feature  in  control  of  structural  acoustics. 


4.2.  Modal  controllability  and  smart  structures  [J5,  C7] 

The  response  of  the  controlled  structure  (see  Fig.  9)  is  governed  by: 


t  t 

Mx  +  Dx  +  Kx  +  jRit-  x)x(x)dx  =  jB^it-  x)u(x)dx 
0  0 

Actuator  ^ 
locations  Control 

&  dynamics 


Structural 

dynamics 


Smart  material 
and  constrained 
damping  layers 


A  condition  on  controllability 
using  modal  information  is 
obtained  as  follows: 

To  control  the  kth  mode, 

where  Ba(s)  is  the  Laplace 
transform  of  the  operator  Bg(t), 
and  and  are  the  eigenvalue 
and  corresponding  mode  shape 
that  are  determined  by  DTFM. 


Active  Constrained 


u(t) 
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Uniqueness  of  the  Modal  Condition 

•  Explicit  relationship  between  controllability  and  modes  to  be  controlled,  and  quantitative 
measure  of  controllability  by  ||v*|| 

•  Applicable  to  structural  systems  with  uncertain/unidentified  parameters  (e.g.,  inherent 
damping,  coupling  of  multiple  bodies) 


Utility  in  Smart  Structure  Design 

•  Optimal  placement  of  actuators  and  active/passive  layers,  and  optimal  shape  design  of  smart 
material  and  damping  layers 


•  Systematic  design  of  feedback  gain  pG(j)  based  on  the  sensitivity  coefficient  of  the  close-loop 
poles  pj, 


dpi^  _  Change  in  closed  -  loop  pole 
d\l  Change  in  control  gain 


—  Q^k  ^  ^(^k  O^k  )^A: 


where  the  operator  Bj(s)  describes  the  locations  and  dynamics  of  the  sensors,  and  G(s)  is  a 
controller  transfer  function  to  be  designed. 


4.3.  Tuned  distributed  vibration  damper 


A  novel  concept  of  tuned  distributed  vibration  damper  (TDVD)  is  introduced.  A  TDVD  consists 
of  an  elastic  continuum  with  distributed  inertia  and  stiffness,  and  a  boundary  conditioner  with 
tunable  parameters.  The  boundary  conditioner  can  be  of  various  forms,  including  conventional 
viscous  damper,  active  feedback  controller,  combination  of  passive  and  active  components,  and 
smart  materials.  The  main  function  of  the  boundary  conditioner  is  to  deliver  passive/active 


damping  for  energy 
dissipation,  while  the  main 
function  of  the  bar  is  to 
amplify  or  augment  the 
damping  delivered  by  the 
boundary  conditioner.  The 
preliminary  analysis  shows 
that  the  TDVD  can  increase 
the  damping  in  a  controlled 
structure  by  several  orders 
of  magnitude,  in  a  wide 
frequency  rggion;  see  the 
example  in  Fig.  10. 

Continued  study  of  this 
interesting  topic  may  lead  to 
new  ways  for  intelligent 
control  of  structures  and 
equipment  in  various 
military  applications. 


Fig.  10  Frequency  responses  of  a  clamped-clamped  beam 
controlled  by  a  TDVD  whose  boundary 
conditioner  is  a  damper,  and  by  a  pure  damper 
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